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The Gaussian chain model is the classical description of a polymeric chain, which provides the 
analytical results regarding end-to-end distance, the distribution of segments around the mass cen¬ 
ter of a chain, coarse grained interactions between two chains and effective interactions in binary 
mixtures. This hierarchy of results can be calculated thanks to the alpha stability of the Gaussian 
distribution. In this paper we show that it is possible to generalize the model of Gaussian chain 
to the entire class of alpha stable distributions, obtaining the analogous hierarchy of results ex¬ 
pressed by the analytical closed-form formulas in the Fourier space. This allows us to establish the 
alpha-stable chain model. We begin with reviewing the applications of Levy flights in the context of 
polymer sciences, which include: chains with heavy-tailed distributions of persistence length, poly¬ 
mers adsorbed to the surface and the chains driven by a noise with power-law spatial correlations. 

Further, we derive the distribution of segments around the mass center of the alpha-stable chain 
and the coarse-grained interaction potential between two chains is constructed. These results are 
employed to discuss the model of binary mixture consisting of the alpha-stable chains. On what 
follows, we establish the spinodal decomposition condition generalized to the particles described by 
the shape of alpha-stable distributions. This condition is finally applied to analyze the on-surface 
phase separation of adsorbed polymers, which are known to be described with heavy tailed statistics. 

PACS numbers: 36.20.-r, 05.40.Fb, 68.35.Dv, 82.35.Gh 


I. INTRODUCTION 

While the theory of Flory provides an accurate de¬ 
scription of the ideal polymeric chains [1], factors such 
as complex environment interactions, adsorption or de¬ 
signed chemical composition can lead to significant devi¬ 
ations from this model. The Flory approach is based on 
the Gaussian chain model, in which the conformation of 
a chain is equivalent to the trajectory of a particle under¬ 
going the thermal Brownian motion [T] . In this model the 
chain is characterized by the Gaussian distribution of the 
nearest neighbor distances, a fact that leads to the entire 
hierarchy of analytical results. In particular, the Gaus¬ 
sian shape propagates to such characteristics as: end- 
to-end distance distribution [1], distribution of segments 
around the mass center of the chain [5] and the coarse¬ 
grained interaction potential between two chains in terms 
of the distance between their mass centers [3]. Deriving 
all of these characteristics is possible due to a single fact: 
the Gaussian distribution is alpha-stable. Since there ex¬ 
ist the entire class of alpha-stable, heavy-tailed distribu¬ 
tions [4] , this suggests that a natural and equally prolific 
generalization of the Gaussian chain model can be based 
on the alpha-stable distributions. Indeed, in this paper 
we discuss the alpha-stable chain model and calculate all 
of the characteristics analogous to the Gaussian model. 

The first goal of this paper is to establish the physical 
context in which the alpha stable distributions are rele¬ 
vant for the polymer sciences. Since the application of 
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alpha-stable distributions (or Levy walks/flights) in this 
context is not an entirely new concept, in Sectionjl^we re¬ 
view the relevant literature. In addition, we provide our 
own simulations of a model polymeric chain under the 
spatially correlated noise, which establish another con¬ 
text for our considerations. 

The main part of this paper is focused on deriving 
and analyzing the different aspects of the alpha-stable 
chain model. In Section mi we introduce the model it¬ 
self. In Section ||V] the distribution of nodes around the 
mass center of a chain is calculated and in Section fVl the 
coarse-grained model of interaction between two chains is 
established. All of these results are analytical and closed- 
form in the Fourier space. 

Another goal of this paper is to analyze the stability 
of binary mixtures composed of the alpha-stable chains. 
Understanding the behavior of binary mixtures is a vital 
problem in industry, medicine, wet nano-technology and 
biophysics. The stability of solution is governed by effec¬ 
tive interactions [5], predicting which is a classical prob¬ 
lem of soft matter physics [5]. In the context of binary 
mixtures of Gaussian particles, Bolhuis et al. found via 
simulations that the interaction between particles of one 
species has also a Gaussian profile, but with an addition 
of a shallow attractive tail [7]. Similar results were pre¬ 
dicted half-analytically via closure-relations techniques in 
UM- On the other hand, a simpler, but entirely analyt¬ 
ical method has been recently proposed by the present 
authors, in HO]. Therein, we have studied the stabil¬ 
ity of Gaussian particles mixtures and our results proved 
similar to the spinodal decomposition analysis [Hill], 
namely, there exist a well defined region of mixing and 
demixing, dependent on the proportion of gyration radii 
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describing the two types of particles. Since our method¬ 
ology from [10] can be conveniently extended to the par¬ 
ticles described with the alpha-stable profiles, we apply 
this approach in the current paper. As the main result 
of Section VI we generalize the spinodal decomposition 
condition for Gaussian particles to the entire class of par¬ 
ticles based on the alpha-stable distributions. We discuss 
the validity of our methodology in Section [vni 

Finally, in Section [VIII| we employ the results regard¬ 
ing phase separation in binary mixtures to analyze the 
phase separation of adsorbed polymers versus their be¬ 
havior in the bulk. As a result we predict the parame¬ 
ters for which the different combinations of simultaneous 
mixing/demixing on the surface and in the bulk can be 
achieved. 


II. LEVY FLIGHTS IN POLYMER SCIENCES 

Except for the Gaussian case, the asymptotic behavior 
of the alpha stable distributions is of the power-law type 
[3] . A random walk characterized by such a heavy-tailed 
distribution of steps is known as Levy flight. It is usually 
difficult to interpret Levy flights in physical terms, thus 
let us discuss three situations justifying such statistics in 
the context of polymers. 

The first scenario can be related to the non-Gaussian 
distribution of segment persistence lengths. The Gaus¬ 
sian chain model is usually interpreted as a continuous 
limit of a discrete model, in which all segments have 
the same persistence length [T]. However, it can be also 
seen as the model for chain made of unequal segments, 
whose persistence lengths follow the Gaussian distribu¬ 
tion. This can be further generalized to the alpha-stable 
distribution, the idea suggested by Moon and Nakanishi 
in [13] . They proposed the Levy walk chain model, based 
on the formalism of turbulent flows m and calculated 
Flory exponents for this model. While no direct exper¬ 
imental confirmation of this idea is known to authors, 
non-Gaussian persistence length distribution might be 
the result of the varying chemical composition of a chain. 
E. g. DNA double strand is characterized by the persis¬ 
tence length much greater than a single base pair m, 
but also certain sequences of chemical monomers can as¬ 
semble into relatively long and stable structures of sig¬ 
nificant persistence length, such as protein domains m- 
A possible realization of such ’Levy flights’ could be the 
intrinsically disordered proteins, in which second-order 
structural motifs such as alpha-helices coexist with dis¬ 
ordered loops [muH]. However, it should be mentioned 
that some numerical experiments on the structure of par¬ 
tially unfolded proteins indicate that Gaussian statistics 
is rather robust m- 

Another scenario is similar to the problem of a 
tracer which mixes one- and three-dimensional diffusion. 
Such motion has been observed experimentally in DNA- 
binding proteins [20] and its simulations revealed the 
heavy-tailed distributions of steps along the polymer in 


certain configurations m- This behavior can be effi¬ 
ciently modeled with Levy flights [22]. In the context of 
polymers, we consider the adsorption of a chain to the 
surface. This problem was first analyzed by de Gennes 
from the scaling perspective in (231123|. Eor the inter¬ 
mediate attraction strength, only some fraction of seg¬ 
ments is attached to the surface, while the loops that 
connect those segments diffuse into the bulk. Consid¬ 
ering the projection of the chain on the surface, it has 
been argued by Bouchaud and Daoud [25| that the pla¬ 
nar trajectory connecting adsorbed nodes can be mod¬ 
eled as Levy flights, since the subsequent adsorbed seg¬ 
ments connected by a loop can be found at abnormally 
long distances. Bouchaud and Daoud predict that for an 
adsorbed polymer its gyration radius parallel to the sur¬ 
face scales as i?g_|| cx for the strong adsorption and 
cx for the weak adsorption |2S|. Within the 

Elory-type theory for Levy flights discussed in [25], the 
former translates into the characteristic exponent of the 
distribution a = 1 and the latter demands a = 2. This 
means that while in the strong adsorption regime the 
Levy flights and the Gaussian statistics are equivalent, 
for the weak adsorption limit should be modeled 

with power-law distributions. 

One final interpretation can be related to the situation 
in which a polymer experiences a random, though spa¬ 
tially correlated behavior of surrounding environment. 
Such conditions occur in the glassy state, in which the 
correlations are exponential [26fl28| . or near crystalliza¬ 
tion, in which case the scale-free behavior results in the 
power law correlations [22]. In [SD] we have simulated 
a two dimensional model polymeric chain driven by the 
spatially correlated noise and observed the effect of spon¬ 
taneous chain unfolding, i.e. a significant number of seg¬ 
ments tends to form linearized structures, scattered along 
the chain. As we have shown in EH, this effect was 
mainly due to the short (2-3 segments) structures, but 
structures up to 50 segments were also observed. Such 
elongated fragments may act as the Levy flights, provided 
that their distribution is wide enough. 

The model from |30| consisted of a bead-spring chain 
with a global Lennard-Jones potential assigned to each 
bead and a second-nearest neighbor harmonic interac¬ 
tion to induce non-linear conformations. The system was 
driven by the noise which spatial correlation function 
read: 


< |(r)4(r-b Ar) >cx exp(—|Ar|/A) (I) 

Eor the purpose of the current paper we have repeated 
the simulations from [30] . replacing the exponential cor¬ 
relations in the noise with a heavy-tailed function, based 
on the Gauchy distribution, namely: 

< KOKr + Ar) >cx 1/ (l -b (|Ar|/A)2) (2) 

The data regarding the linearized fragments has been 
gather in the same fashion as in m- To improve statis¬ 
tics, the single set of parameters was simulated 128 times. 
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otherwise, the details of the simulations remained the 
same as in [30] and |3T]. In Fig. we include the rep¬ 
resentative probability distributions S'„ of finding the n- 
segments long structure in the chain geometry. The data 
has been gathered from the regime of noise dominated 
dynamics. With the growing correlation length A the 
distribution Sn gradually develops a linear region in the 
log-log plot. For A > 40, where relevant, we have fitted 
Sn for n > 17 with a power-law model 
We obtain a ranging from 1.18±0.39 to 2.44±0.54, with 
the relative error typically at the level of 20-30%. The 
uncertainty intervals for these values of a overlap with 
the interval a G [0,2], which is expected for the asymp¬ 
totic behavior of the alpha-stable distributions |4]. For 
comparison, the data have been also fitted with the expo¬ 
nential decay model Sn = . While both the linear 

and exponential fit describe the tail part of similarly 
well (in both cases B? ~ 0.8) and, most probably, even 
for A = 50 the distribution Sn eventually develops the 
exponential decay for n ^ 50, a ’power-law-like’ region 
suggests that in the special conditions of long-range spa¬ 
tial correlations the alpha stable distributions might be a 
more relevant description of the chain statistics than the 
Gaussian distribution. For comparison we also include in 
Fig. [^the data from |30j , which preserve the exponential 


form in the entire range of parameters. 


III. ALPHA-STABLE CHAIN MODEL 


In the Gaussian chain model, the geometry of a chain 
is described as a random walk trajectory, in which the 
distribution of the distances between nearest neighbors 
is Gaussian, namely [Tj: 


G(|r,+i -r,|) = 


f2TTb‘^\ f £)(ri+i-r,)2 

exp -■ 


V D 


262 


( 3 ) 

Here, D is the dimension of the system, 6 is usually in¬ 
terpreted as the length of a segment and is the vector 
position of i-th node. The characteristic function of G 
reads: 


(/>G(fc) = exp (4) 

Let us now consider the generalization of G(|ri+i — rij) 
to the alpha-stable distribution P(|r.i+i — r.ij). The mul¬ 
tivariate alpha-stable distributions are defined in terms 
of their characteristic functions, which can be written in 
the following parametrization 


exp - 


0(k) = 


j |k • s|“ ^1 — isgn(k • s) tan F((is) -t- ik • /X j for a^l 
^ Sd / 


exp 


fSn 


k • s| ( 1 -f JSgn(k • s) 


2 In k • s 


( 5 ) 


F(ds) + ik ■ fj,\ for 0=1 


In the above definition r((is) stands for the spectral mea¬ 
sure defined on the D-dimensional unit sphere Sd,^-s 
denotes the scalar product and /x is the vector mean 
value. Since we are interested in the spherically symmet¬ 
ric distributions, we choose the uniform spectral measure 
F(ds) = const [35]. Under such choice, and assuming 
/X = 0, the general parametrization can be simplified to 
the following form: 


(p(k) = exp {—ck°‘) 


( 6 ) 


where fc = jkj and c is a constant. Eq. Q is a special 
case of (j){k) and our choice of c should agree with (j)G{k) 
for a = 2. Therefore, we postulate that c reads: 


26“ 

1d~ 


( 7 ) 


and, finally, the nearest neighbor spatial distribution in 
the alpha-stable chain model reads: 


P(|ri+i-ri|) = 


1 

(27r) 


dkexp ( ik • (r^+i — r^) — 

( 8 ) 


Having established the nearest neighbors spatial distri¬ 
bution P(|ri+i — Til), we can calculate such distribution 
for any pair of segments, namely: 


-Pdu-rjl) = 

= J dri+i... 
1 


(27 




/ dvj-i P(|r„ - r„_i|) = 

J dk exp ^ik • (ri - r^) - ^1* - 


( 9 ) 


where we made use of alpha-stability. For a = 2 this for¬ 
mula comes down to the well known result for Gaussian 
chain [1]: 


G(|ri-rj|) = 


2n\i-j\b'^\ 


D 


exp 


D{ri - Vjf 
2|* - j|62 


( 10 ) 
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Figure 1. The distribution of linearized polymer structures 5„ in the system driven by the spatially correlated noise. Each 
plot represents the results for noise amplitude 2kT and correlation length A. Triangles: exponential noise correlation function 
< 4(i’)$(r + Ar) >oc exp(—|Ar|/A), dots: Cauchy noise correlation function < ^(r)^(r + Ar) >oc 1/ (l + (|Ar|/A)^). The tail 
behavior (n > 17) has been fitted with exponential decay model Sn = (dashed lines) with bi given on the plot. For 

A > 40 the data has been also fitted with power-law S„ = for n > 17, indicating the power law asymptotic behavior. 


It is also possible to calculate a = l case explicitly: 


Pa=i,D=3{\Yi - r^l) 


167r|z - j\ 
3 


1 


Mi-j? 



( 11 ) 


In particular, taking as i and j the first and the last 
segment respectively, we obtain the end-to-end distance 
distribution. 

A classical problem in polymer physics is to predict the 
scaling behavior of radius of gyration Rg with growing N. 
In order to do so, we will calculate Rg using a method 
mentioned in na, namely Rg =< r“ where < . > 

denotes the average: 


pletely factors the dependence on b and N out of the 
integral. Therefore, the scaling reads Rg (x 


IV. DISTRIBUTION OF SEGMENTS AROUND 
THE CENTER OF MASS 

The central problem of deriving the coarse-grained in¬ 
teraction between two chains is to find the distribution 
of segments around the mass center of a single chain. Let 
us consider an V-segments long chain, described by the 
nearest neighbor probability given by (|^. The position 
of the mass center reads: 


Rg — 



( 12 ) 


This result is obtained via the change of variables k' = 
6(2V/Z7)^/“k and r' = b~^(2N /which com¬ 


1 ^ 

(13) 

The probability that each segment occupies its position 
i/f*i under condition that the mass center is positioned at 
R reads: 

iv-i / 1 ^ \ 

P(ri,..., r^vIR) = n ^(l’'ui - 1‘il)'^ II - H 

(14) 
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From this expression we can calculate the probability of 
finding j-th segment at some position relative to the mass 
center: 

P(|rj-R|)= J dri...J dr^_i J drj+i...J drArX 


N-l / N 


N 


(15) 


The integrals in the above expression can be done par¬ 
ticularly easily, if we switch to relative variables: 


Ari_j =r,- ri_i for i > j 
Ari_j = r^_i - Ti for i < j 


(16) 


which allows us to express as: 

N-i 


Ar+„ for N > i > j 

n—1 

j-'i 

y^Ar -n for 1 <i < j 


(17) 


n—1 


In these coordinates the position of mass center reads: 

i—1 n—1 n—1 

^ ( 18 ) 

The change of variables (161 is linear and its Jacobian 


is equal 1, so applying the new coordinates to (15), we 
obtain: 

N-j . 

P{\rj - R-l) = n / dArr,P{Arn)x 


n—— 7 
n/0 

N-j 


Further, we integrate out every component of Ar±„, 
which introduces multiple Dirac-delta functions: 

Pi\r,-K\) = 


(27r) 


D ' 


/N-j 


= I I rfk+n'('(fc-Sn)<5 ( k+n “ ^ ^ ^ ^' ^0 


^ n—1 


0-1 


/ dk_ricj){k-„)S ( k_„-—k| 


\n—l ‘ 


N 


•^-^ko ] 1 


/N-j 


(27r)-® 

x(n< 

\n=l 


dko I 

\n=l 

j -n 
N 


0 e 


N — j — n + 1, 
- - - 

N 

iko'(R-rj) 


( 21 ) 


At this point we apply the explicit form of (f>{k), so the 
final expression for P{\rj — R|) reads: 


Pi\rj - R|) = 


1 


(27r) 


D 


dk-o exp 


26“ fN-j-n+1 


D 


N 


iko • (R- r^)- 

i-i 


E 

n—1 


J - n 
N 


^0 

( 22 ) 


Expression (22) gives the probability of finding j-th seg¬ 
ment in the vicinity of mass center, so the probability of 
finding any segment reads: 


PcM(|r-R|) = ^y]P(|r,-R|) (23) 


N 


N 


i=i 


X (5 R - -f y]] 


n—1 


N — j — n + 1 
N 


J - n 


+ E 


n—1 


where the factor 1/A provides normalization. Let us as¬ 
sume now that A is a large number, so both n/A = q and 
j/N = q' can be treated as continuous variables, hence 
we can simplify: 




■Ar_„j 

N-j 

E 

(19) 

n—1 


A 


The following step is to express P(Ar„) in (19) in terms 
of its characteristic function : 


1 / r f 

-R|) = / 27 r)^-P n / J c^k„e*‘""'^'’"(j(fc„)x 

' ' n=—j 

n/0 

X J dko exp ^iko • (R — Yj)+ 

AT ■ ,1 


Y(1-9') 




(24) 


i-i 

E 

n—1 


J - n 
A 


/■9 N 

A dq{q' - qr = (25) 


a -I- 1 


The final expression for the distribution of any segment 
around the mass center reads: 


PcM{\r -R\) = J dq' J dkox 


E N — j — 77,-|-l . J 

---ko • Ar+n + ^ ^ ^ I'-r 


J-n. 


n—1 


n—1 


A 


X exp 


(20) 


9 

*ko • (r - R) - ° ((1 - 9')“+' + 9'“+') 

D{a + 1) 


(26) 

































6 


While this expression is exact, the integral with respect 
to q' makes it unwieldy. We can simplify it by resorting 
to the integral mean value theorem. Namely, there exists 
such Qq G [0,1], that: 


PcM(|r-R|) = 
exp jko•(r — R) — 


1 


(27r)^ 

2Nb°‘k 


dknx 


D{a + 1) 




(27) 


On the other hand, since the integrand of (26) is a peak 
function of q' with maximum at q' = 1/2, this value 
contributes the most to the integral. For this reason, we 
will approximate ( |27[ ) by imposing Qq = 1/2. Once again, 
this results in the characteristic function of Pcm (|r — R|) 
in a form exp(—c(a)fc“), where: 


c{a) 


2 N ny 

D{a + 1) UJ 


(28) 


V. COARSE GRAINED INTERACTION 
BETWEEN TWO NON-GAUSSIAN POLYMERS 

Having established PcM(|r — R|) for a single chain, 
we can calculate the coarse-grained interaction between 
two chains in terms of the distance between their mass 
centers. From now on, lower index numerates the type of 
particle, so PcM.i (Ir—Ril) describes an i-th type of chain 
characterized by Ni segments, constant and exponent 
ai- We can follow the reasoning of Flory [3] and assume 
that the systems suffers an energetic penalty if a seg¬ 
ment belonging to the first chain invades a small volume 
in the vicinity of a segment belonging to the other chain. 
For a single site r, the probability of such event is propor¬ 
tional to PcM.i(|r — Ri|)PcMj(|r — RjDdr. Therefore, 
the entire interaction reads: 


that c = Ci{2) + Cj{2) and hence Cij = (47rc)^fo. This 
can be generalized for ai = aj = a by: 

Cij = (47r(c^(a) -k Cj(a)))'^/“ (31) 

For the case of ai ^ aj the constant Cij cannot be 
uniquely deduced from the dimensional analysis, thus we 
will restrict our further considerations to the potentials 
with common a. 


VI. EFFEGTIVE INTERACTIONS AND 
MIXTURE STABILITY 


Once Vij has been found, we can analyze the inter¬ 
actions in binary mixtures. The system is described by 
three microscopic potentials in the form (29), where Vu 


and V 22 are the internal interactions of each species and 
V 12 is the cross-species interaction. When the behavior of 
one species in a mixture is considered, the presence of the 
other species modifies the microscopic interaction mm- 
The additional potential, known as the effective interac¬ 
tion, is of entropic origin [SJIS] and it is a key-factor in 
controlling mixture stability and demixing. The predic¬ 
tion of effective interactions from arbitrary microscopic 
potentials is usually a challenging numerical task, but in 
our previous work uni we have proposed a simple ana¬ 
lytical method, suitable for soft interactions. According 
to [To], the effective interaction can be estimated by: 


UeffiAR) = - 


I 


(27r) 


D 


dke 


ik-AR 


|Vl2(fc)l^ 

V22{k) 


(32) 


where Vij(k) = f c?rexp(zk • r)Vij{r) is a Fourier trans¬ 
form of relevant Vij (r) and H is the volume in the recipro¬ 
cal space. Substituting (29) with relevant constants given 


by (311 into the expression for effective interactions, one 


obtains: 


U„(|R,-R,|) = 

= CijCij j drPcM,i(\^-'^i\)PcM,ji\r-Rj\) = 

= J dkexp (ik- (R^ - Rj) - Ci(ai)k°‘' - Cj{aj)k°'^) 

(29) 

Assuming that has a dimension of energy, it is neces¬ 
sary to introduce an additional constant Cij , which has a 
dimension of volume. We can deduce this constant from 
the case of = aj = 2, for which we obtain the fol¬ 
lowing ’universal’ Gaussian potential |5] and its Fourier 
transform: 

U(r) = eexp V(k) = e (Attc)^^^ (30) 


UeffiAR) 


1 ei2 (47r)^/“(ci(Q;) -|- 02 ( 0 ;))^'°^ 

" (^^ (2c2(a))^/“ "" 

Jn 

(33) 


The total interaction for the first species in the mixture 
reads: 


UtotiAR) = Un(Ai?) + UeffiAR) (34) 

or explicitly: 

Utotir) = 

( ,, ( iiD/a e ?2 (47r)'°/“(ci(Q:)-kC 2 (a))^^/“\ 


When e is independent from V, this potential is perceived 
as an accurate and reliable model for interaction of iden¬ 
tical chains iniiTj. Comparing (j^ to (|^, one can see 


(27r) 


D 




(35) 
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One can see that Utot{^R) and Vii(A_R) have the same 
shape, up to the scaling constant S: 


S = eii (47rci(a))^/“ 


ei2 (47r)^/“(ci(a) + C2{a)f‘^/' 
£22 (2c2(a))'°/“ 


(36) 

S has a complicated form and it can take both the nega¬ 
tive and positive values, depending on the parameters. 
The change in the sign of the total interaction indi¬ 
cates a remarkable change in the behavior of the system. 
Namely, for S' > 0 the interaction between the particles 
of the first species are purely repulsive, which means that 
these particles will disperse in the volume. Conversely, 
for S < 0, the first species of particles interacts via at¬ 
tractive potential, which results in the clustering of these 
particles and demixing in the system. Equating S to 0, 
the condition for demixing reads: 


£11622 ^ / (ci(q) -I- C2(a))^ \^^° 
e|i V 4ci ( 0 ) 02 ( 0 ) ) 


(37) 


Let us now analyze the condition (371 and introduce a 
common energy scale: 


e = 


£12 


and: 


g = 


-v/£11£22 

1 1 \ l/“ 
ci(a) ^ 

C 2 (a) 


N 2 


l/a 


for which condition 


£ > 


can be reduced to: 


4<7° 


(l + g “)2 


D/(2a) 


(38) 


(39) 


(40) 


Recalling the equation (12) for the radius of gyration Rg, 


one can see that for chains characterized by the distribu¬ 
tions sharing the same o, the parameter g becomes the 
ratio of Ra'. 


9 = 


R, 


9.2 


(41) 


For 0 = 2 and D = 3 the condition (37) becomes ex¬ 


actly the spinodal decomposition condition for Gaussian 
particles, as given in m and m, namely: 


£ > 


25 


l + (?2 


3/2 


(42) 


Therefore, the condition (40) is a direct generalization 


of the spinodal decomposition to the systems of particles 
described with alpha-stable distributions. 

The condition (40) is plotted in Fig. ^ For every pair 
of g and a its value fits between 0 and 1. In the entire 
range of a, the region of mixing (below the surface) pre¬ 
serves the features of the Gaussian case, namely it falls 



Figure 2. The spinodal decomposition condition ( |40[ ) for the 
binary systems of particles described with the alpha-stable 
statistics as a function of gyration radii ratio g = Rg,i/Rg ,2 
and the distribution exponent a for D — 3 dimensional sys¬ 
tem. e = ei 2 /-^eiie 22 is the common energy scale. In the 
region above the plotted surface the binary system undergoes 
demixing due to the prevalence of the attractive effective in¬ 
teractions, in the region below the surface the effective inter¬ 
action is repulsive and the system is homogeneous. 


rapidly to 0 for <C 1 , reaches the single maximum at 
g = 1 and asymptotically decreases to 0 for g I. How¬ 
ever, as a decreases to 0 the mixing region for 5 I be¬ 
comes wider, asymptotically reaching the region defined 
by £ > I. This means that the gyration radius ratio g 
becomes less and less relevant for the mixing of chains 
characterized by very wide distributions. The changes in 
the mixing region shape are much more pronounced for 
a < 1. 


VII. DISCUSSION 


It is known that the spinodal decomposition condi¬ 
tion for Gaussian particles leads to the predictions on 
mixture separation which are qualitatively and quanti¬ 
tatively comparable to the more advanced methods )I 2 j . 
Thus, a similar efficiency can be expected from (40), at 


least for a mildly deviating from 2. However, some pos¬ 
sible issues should be commented. 


The fact that the total interaction (34) can become 


entirely negative is unrealistic. This evident problem is 
mitigated by the fact, that the energy density of a pair 
interaction behaves as Utot{r)r^~^dr. Therefore, 
factor suppresses the lack of repulsive core at short dis¬ 
tances and amplifies the influence of the attractive tail. 
The unrealistic shape of Utot is also a consequence of 
the way the potential Vij{r) given by (29) is constructed. 
This potential is mean-field in its nature and its width is 
governed by the constant Ci{a) + Cj{a). In the context 
of Gaussian particles, while the dominant shape of the 
interaction between two separate chains is agreed to be 
























Gaussian [71 [53], there is an open problem of whether 
there are additional components |38| or how the width 
of such Gaussian is related to the gyration radii of the 
component chains m- A similar problem is relevant in 
our case and the choice of the width constant different 
from Ci(a) + Cj(a) might result in a more realistic shape 

of [ftot 


VIII. PHASE SEPARATION IN THE 
ADSORPTION OF GAUSSIAN PARTICLES TO 
THE SURFACE 


The result (401 is particularly interesting in the context 
of the already mentioned work of Bouchaud and Daoud 
[25], where the gyration radius parallel to the surface is 
calculated for an adsorbed polymer. As mentioned in 
Section [IT] the characteristic exponent for the distribu¬ 
tion on the surface reads ass = 2 in the strong adsorp¬ 
tion limit and a^s = 1 in the weak adsorption limit |25| . 
Considering the adsorption from the binary mixture, (40) 
provides the condition for homogeneous versus inhomoge¬ 
neous adsorption, i. e. the ratio of gyration radii parallel 
to the surface (.Rg,||) decides whether both species cover 
the surface in a homogeneous manner or they separate 
into the islands consisting of the particles of the same 
type. On the other hand (40) allows us to compare for 
which parameters the separation on the surface and in 
the bulk co-appear. 

Let us consider a simple binary system in which the 
behavior of chains in the bulk (Uf, = 3, at = 2) is Gaus¬ 
sian, but on the surface it is characterized by Dg — 2 
and Oss or a^s- The types of particles differ by the num¬ 
ber of monomers Ni and their persistence length bi. The 
condition (37) reads: 



In the strong adsorption limit it is always true that 
fss,s > fb,b for any 61/62 and N1/N2. Therefore, as¬ 
suming that in this system e is the same on the surface 
and in the bulk, three scenarios of mixing/demixing are 
allowed. First, for (i) „ > f^ i, > e the solution in the 

bulk is homogeneous and so is the coverage of the sur¬ 
face. Conversely, for (ii) e > fss,s > fb,b the separation 
is simultaneous on the surface and in the bulk. Finally, 
for (iii) fss,s > e > fb,b demixing in the bulk occurs, but 
the surface coverage is still homogeneous. 

In the weak adsorption limit {a^s = 1), the situa¬ 
tion is more complicated because both fujs,s > fb,b and 
fws,s < fb,b are possible, depending on 61 / 62 , Ni/N 2 . 
Replacing fss,s by fws,s in the inequalities from the pre¬ 
vious paragraph one obtains the conditions for separation 
scenarios (i)-(iii) in the weak adsorption limit. However, 
there exists the additional region in which it is possible 
that (iv) fb^b > e > fws,s- In this region the separation on 


Phase separation condition ^ and fb_b for di/b2=0.1 
€ 



bulk fbb- Db=3, 0^=2 


Figure 3. The exemplary comparison of the phase separation 
conditions on the surface in the weak adsorption limit (/ms.s, 
solid line) and in the bulk {fb,b, dashed line) for the binary 
system of Gaussian polymers characterized by the persistence 
length ratio 61/62 = 0.1. f-ws,s and fb,b are defined by (431. 
fws,s and fb,b divide the plot into four regions: (i) simultane¬ 
ous mixing on the surface and in the bulk, (ii) simultaneous 
demixing on the surface and in the bulk, (iii) mixing on the 
surface, demixing in the bulk, (iv) demixing on the surface, 
mixing in the bulk. 


the surface occurs, while the solution in the bulk is still 
homogeneous. In Fig. [^the exemplary phase separation 
diagram for the weak adsorption limit and 61/62 = 0.1 
is presented, which contains all of the phase separation 
scenarios (i)-(iv). Fig. [^shows the difference fws,s — fb,b, 
which indicates where scenarios (iii) and (iv) are allowed. 
In particular, for 61/62 —t 0 and 61/62 1 the scenario 

(iv) becomes almost inaccessible, while it is allowed in 
the vicinity of 61/62 — 1 . 

While these considerations show that the behavior of 
mixture can be designed by the choice of 61/62 (which 
is dependent on the chemical composition) and N 1 /N 2 , 
some additional effects might affect the result. One is¬ 
sue is that there exists the additional entropic attraction 
between surface and the bigger particles 0131] which is 
not included here. Another problem is that our consid¬ 
erations are relevant for the thermodynamic equilibrium, 
which might be difficult to achieve in the adsorption pro¬ 
cess. 


IX. SUMMARY 

In summary, in this paper we have presented the gen¬ 
eralization of the results known from the Gaussian chain 
theory to the particles described with the alpha-stable 
distributions. As expected, it is possible to obtain a sim¬ 
ilar hierarchy of analytical results ranging from end-to- 
end distribution up to the effective interactions in binary 
mixtures. Typically for alpha-stable distributions, we ob- 
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tained the closed-form formulas in the Fourier space. Our 
theory also allows us to generalize the spinodal decom¬ 
position condition from Gaussian particles to the alpha- 
stable particles. This can be readily applied to the prob¬ 
lem of mixing/demixing of adsorbed polymers, as we also 
show. Our results might be further utilized in the context 
of Levy flights applications reviewed in Section [Tl] 


Figure 4. The difference fwa,a — fb,b indicates which sur¬ 
face/bulk demixing scenarios are allowed (see Section EM 
and Fig. [^for explanation). White meshed region: /■ws ,3 — 
fb,a < 0 indicates scenario (iii) allowed (demixing on the sur¬ 
face, mixing in the bulk), complement region: f-wa,a-fb,b > 0 
scenario (iv) allowed (mixing on the surface, demixing in the 
bulk). 
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